New mechanism of pulsar radio emission 
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It is shown that pulsar radio emission can be generated effectively through a streaming motion in the polar-cap 
regions of a pulsar magnetosphere causing nonresonant growth of waves that can escape directly. As in other 
beam models, a relatively low-energy high-density beam is required. The instability generates quasi-transverse 
waves in a beam mode at frequencies that can be well below the resonant frequency. As the waves propagate 
outward growth continues until the height at which the wave frequency is equal to the resonant frequency. 
Beyond this point the waves escape in a natural plasma mode (L-0 mode). This one-step mechanism is much 
more efficient than previously widely considered multi-step mechanisms. 

PACS numbers: 97.60.Gb,97.60.-s 

Since the discovery of pulsars in 1967 the mechanism of generation of their highly nonthermal (brightness temperatures up 
to 10 29 K) pulsed radio emission (in the range 10 s — 10 11 Hz) remains one of the most intriguing astrophysical puzzles [jj]]. 
Compact sizes (radius R p ~ 10 6 cm), fast rotation (period P ~ Is), and superstrong polar magnetic fields (B ~ 10 12 G) result 
in the efficient avalanche production of an ultrarelativistic pair plasma (Lorentz factor 7 p ^> 1) in the vicinity of the magnetic 
poles of the neutron star (see, e.g. Ref. fy. This plasma flows outward along the open magnetic field lines and escapes the pulsar 
magnetosphere as a relativistic wind beyond the light cylinder, R = cP/2ir. In a standard polar-cap model, a primary beam 
(Lorentz factor 7& 3> j p ) of particles of one species propagates through a secondary pair plasma [JTJ . The superstrong magnetic 
field of the pulsar implies a very short lifetime for the electrons and positrons to radiate away all their perpendicular momenta, 
so that the plasma distribution is one dimensional. The properties of this pulsar plasma determine the natural wave modes, and 
the problem is to explain how excitation of these modes occurs and how it produces the observed radiation that escapes from the 
magnetosphere. 

It is probable that the radio spectrum forms in the inner magnetosphere Q], where the infinite magnetic field approximation is 
appropriate. The properties of low-frequency (well below the cyclotron frequency) waves in a one-dimensional, relativistic pair 
plasma have been extensively studied (e.g., Ref. |] and references therein). It has been found [|[ ^[ that for a rather wide class of 
plasma distributions the natural modes are the electromagnetic i-mode, u> — kc, and two mixed (with transverse and longitudinal 
components of the electric field vector) modes, the almost nondispersive Alfven mode, u> — k»VA (here and hereafter subscripts 
|| and _L refer to the direction with respect to the external magnetic field), and the L-0 mode which has a long wavelength cutoff. 
The L-O mode is superluminal, u>/kn > c, near the cutoff but can be subluminal for sufficiently small angles, 9, of propagation 
at high frequencies. The critical features of any plausible mechanism for the radio wave generation can be summarized as 
follows iQ. Only t- waves (which cannot be generated through a beam instability) and waves which eventually appear on the 
L-0 branch can freely leave the pulsar magnetosphere. Beam instabilities are widely favored, and in the pulsar magnetosphere 
are of hydrodynamical type, where the whole beam excites the modes, in contrast to kinetic instabilities which are driven only 
by a group of resonant particles [Q]. 

One of the most widely favored scenarios for the radio emission mechanism involves a resonant instability in which an 
energetic beam causes quasi-longitudinal subluminal waves to grow. This mechanism encounters certain difficulties. The first 
difficulty is that the emission mechanism is indirect: the postulated growing waves need to be converted into quasi-transverse 
waves (in the t or L-O modes) via some secondary (e.g., nonlinear) mechanism before they can produce escaping radiation. A 
direct mechanism in which the escaping waves are generated directly by the instability would be preferable, and the mechanism 
proposed here has this desirable feature. The second difficulty concerns the beam: when the generation is attributed to the 
primary, highly relativistic, 75 ~ 10 7 , the growth is too slow to be effective. Thus, efficient beam instability requires a denser, 
lower energy beam. One suggestion is that such a beam results from nonstationary avalanche pair generation [g[ ||] . Alternatively, 
even in the steady avalanche regime, the high energy tail of the pair plasma distribution [|KJ] transforms into a dense beam due 
to the inverse Compton scattering [[TT|]. The resulting distribution typically consists of a plasma body with j p ~ 10 1 — 10 2 and 
a beam with 7b ~ 10 3 . In the steady case such a beam forms at about 10i? p . Although existence of such a beam seems plausible 
it is still a model dependent assumption. The third difficulty is related to the nature of the wave generation. Most authors have 
concentrated on resonant excitation of parallel (to the magnetic field) propagating modes (Langmuir waves) [j~2|, |T^|. 

where it is widely believed that there is a sharp maximum of the growth rate. Growth for slightly oblique propagation has 
been considered |Q, [l4|] and also found to favor quasi-longitudinal L-O waves. Only subluminal waves can be resonantly 
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excited by a beam. L-O mode becomes subluminal only at frequencies ~ j p x plasma frequency [g]. Growth is possible at the 
resonant frequency, which is just above the frequency where the waves become subluminal. These L-O waves can eventually 
escape directly. However, the frequencies are too high to account for the broadband pulsar emission extending to much lower 
frequencies (so that a conversion mechanism to lower frequencies would be required). Moreover, growth is restricted to a narrow 
range of heights in the magnetosphere where the resonance condition is satisfied, restricting the growth factor to too small a value 
to allow effective growth [^. These constraints led to our relatively pessimist view in |^] of the effectiveness of direct growth of 
escaping quasi-transverse waves. 

In this Letter we show that a nonresonant version of the instability causes a beam mode to grow over a broad range of lower 
frequencies. These waves can grow over a large range of heights. As the propagate outward that ratio of their frequency to 
the resonant frequency decreases, and at which height at which this ratio become unity, they evlove into L-O mode waves, and 
escape, without any secondary conversion process being required. 

We use the following notation: p = mu is the (one dimensional) particle momentum, with m the electron mass, and v = u/7, 
7 = (1 — v 2 )^ 1 / 2 = (1 + u 2 ) 1 / 2 is the particle velocity in units with c = 1. We assume that the distribution function consists 
of a pair (p) and a beam (b) component, /(it) = n p f p (u) + nbfb(u) with J f p du — J fbdu = 1. We adopt the infinite magnetic 
field limit, which is justified in the inner part of the plasma magnetosphere where the beam-plasma system is formed [jilt], and 
then it is not necessary to distinguish between electrons and positrons, which contribute in an identical manner. The dispersion 
relation for this beam-plasma system takes the following form P, pL 0]: 
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with z = w/k\\, fc|| = fccos6>, k± = ksin9, and where e\\ = 1 — (lu 2 / k 2 )W (z) + e& is the parallel dielectric constant. 
Here uj p — (87m p e 2 /m) 1 / 2 (equal electron and positron densities) is the plasma frequency for the pairs, W — (v — z — 
ir)" 1 (dfp/ du)du, r — » +0, is the dispersion function for the pairs; the contribution of the beam is eb = — w 2 7 & ~ 3 (w — k^vt,)~ 2 . 
We consider a cold beam for simplicity; the cold approximation is justified when the instabilities are know to be hydrodynamic. 
The dispersion relation ([[]) takes the form 
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In the absence of the beam the dispersion relation for the natural modes of the pair plasma is K(z) = 0, which can be written in 
the parametric form (recall c = 1) 
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The inclusion of the beam introduces additional solutions, called beam modes, and hydrodynamic instabilities may be attributed 
to a beam mode becoming intrinsically growing. The instability is said to be nonresonant when the beam mode does not coincide 
with a natural mode of the pair plasma, and resonant when it does. The contribution of the beam is significant only when the 
denominator in the left hand side of (||) is small, that is, near z = Vb- Writing to = k\\Vf, + 5u>, \5lo\ <C |w| one finds 
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The unstable solution Sui — iT, T > 0, exists for K(vb) < 0, that is, uj 2 < u> 2 W(vb)v 2 (v 2 — l)/(w 2 — cos 2 9), which requires 
W{v b ) > 0. This is the nonresonant beam instability which sets on the beam mode u> w kuVb. 

The expression (Q) becomes invalid for K(vb) — > 0, that is, where the beam mode uj — k\\Vb resonates with the L-O mode 
(which is the solution of K (vb) — 0). In this case the right hand side of @) should be Taylor expanded up to the first nonzero term 
which immediately gives (see, e.g., Ref. @) T r = Im 8w = (\/3/2)(w^/7^') 1 / 3 , with K' = (dK/dto) res = 2v 2 ^ tan 2 6/w- 
ujpVfrW (vb) I uj 3, , and W'(z) = dW/dz. For a wide class of distributions the approximation W'(vb) ~ 7 2 W / (w; ) ) holds. This 
implies that the resonant growth rate is insensitive to 9 < 7 P /7h, and decreases slowly with 9 > 'Jp/'Yb- This point is important 
in the following discussion: if the growth rate were very sensitive to 9, then a small change in 9 as the waves propagate outward 
along the curved field lines would restrict the distance over which growth can occur, and hence severely limit the possible growth 
factor. 

The polarization of the unstable mode is given by (E±/Eu) = tan#(l — z 2 ) -1 = 7 2 tan(9. Thus, the unstable mode is 
quasi-longitudinal (E || B ) for 9 < r ) b ~ 2 and quasi-transverse (E _L B ) otherwise. The polarization of the growing waves is 
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relatively unimportant in practice: the polarization evolves as the waves propagate outward and the observed polarization may 
be quite different from the polarization at the point of emission [p"5|], 

The derived expressions give the growth rate for arbitrary propagation angle and beam parameters. For small S « 1 and a 
highly relativistic beam OjRil- \j2~f 2 , jb 3> 1, the growth rate for the nonresonant instability (Q) reduces to 



t i -1/2 



(5) 



which is a slightly increasing function of 9. In the range 7^ < 9 < 7^" , where the waves are quasi-transverse, the dependence 
on 9 is negligible. The maximum unstable frequency uj 2 n = uj 2 W '(vb) / (1 + lid 2 ), is also the frequency at which the resonant 
instability occurs, and it is also almost independent of 9 in this range. We conclude that under quite general conditions the beam 
instability generates weakly oblique, quasi-transverse waves with a similar efficiency to that for quasi-longitudinal waves. 

To illustrate the results in the simplest possible way, we use the cold plasma approximation, /„ = 5(u — U p ). In this case one 
has W{z) = r Yp 3 {z — v p )~ 2 . The parametric equations (||) become k 2 = w 2 (z 2 — l)/7p(z 2 — cos -2 9)(z — v p ) 2 , uj — k\\z. For 
small 9 and 7b ^> j p ^> 1, the growth rate of the nonresonant instability becomes 

T„ = 3/2 [4^7 P /- 2 - 1 - ll& 2 ] ~ 1/2 • (6) 

Thus the growth rate for the resonant instability also decreases monotonically with increasing 9. Although these results are 
derived using the cold plasma expressions, as we plan to show in detail elsewhere, they are illustrative of a rather wide class of 
distributions due to the fact that the hydrodynamic instability is insensitive to the details of the beam distribution. 

For small propagation angles 9 < l/ib the growth rates are almost independent of 9, and we can use the following ap- 
proximations. In the ultrarelativistic limit 75 ^> 7 p ^> 1 one has T n = uJblb 3 (ioj^jp/u) 2 — l) -1 / 2 . For low frequencies, 

1/2 3/2 1/2 1/2 

ui -c 2uj p jp , this expression simplifies to r„ = ujujb/2uj p j b j p . The resonant frequency is u> r = 2uj p j p , and the res- 
onant growth rate is T r = 3 1/2 2~ 4/3 (cj p w 2 ) 1/3 /7 b 7p /2 . Note that T n (uj = uj r )/T r — (7 P /7&) 1/2 («&/ n p) 1/3 implies that for 
moderate n p jp/nbib the ratio of the two growth rates is of the order of unity. One can approximate the growth rate in the whole 
range by T — (uj/2^ 2 jp^ 2 )H(uj r — ui), where H(x > 0) = 1 and H(x < 0) = 0. This approximation is also valid if the 
above ratio is small, except in a narrow frequency range around that resonant frequency. We exploit this approximation in our 
estimates below. Note that the condition jb 7p is made for simplicity, and it is not an essential condition for the instability 
to operate. 

The direct excitation of quasi-transverse waves is a fast process, faster than any nonlinear conversion mechanism. Let the 
growth rate be T(uj), which is a function of the plasma parameters, n p , rib, Jp, lb, and n p , through which it depends on the 
radius, R from the center of the star. The wave amplitude evolves according to (da u /dt) = T(ui, R)a u . With the plasma 
streaming outward at close to the speed of light, the solution implies a LU (R) = a u (Ro) exp(/ fl T(uj, R)dR), where i?o is the 
radius where the instability sets in. The power spectrum of the escaping radiation is proportional to the square of this amplitude. 

In a homogeneous static plasma the fastest growing mode is the resonant one, and one might expect that the resonant frequency 
ultimately dominates the spectrum. In the inhomogeneous plasma of the pulsar magnetosphere the conditions change with the 
radius. A wave which is excited at the frequency u> at radius Ro propagates outward into the lower density plasma. A wave 
initially at resonance does not remain resonant as it the propagates. The frequency width of the resonance is ~ r r uj r , 
where r r cx ui r cx n 1 / 2 oc 

j R-3/2_ p or 

a given uj the resonant condition \w — u r (R)\ < T r (R) is satisfied only for a small 
AR/R ~ T r /uj r . As a consequence, the resonant growth condition can be met only for a short time, and effectively only at a 
single height in the magnetosphere, which places a severe restriction on the gain factor G = exp 2 J TdR. It is the gain factor G 
which determines the efficiency of the wave generation and not the local growth rate V. It was for this reason that a pessimistic 
view of the effectiveness of the resonant instability was taken in [Q]. However, a wave at a given uj can grow nonresonantly for 
uj < uj r , and the slightly lower growth rate for the nonresonant instability, compared with the resonant instability, is relatively 
unimportant compared with the much greater distance over which nonresonant growth occurs. The much longer growth path 
through the magnetosphere results in a much larger gain factor. We note that during the propagation the wavevector, in general, 
deviates from the initial propagation direction, so that the propagation angle 9 changes, and in principle this would limit the 
growth if the wave moves out of resonance as 9 increases. This issue will be studied elsewhere. A wave that starts growing 
nonresonantly with w < w r at some radius, Rq, keeps growing while propagating outward until uj r oc R~ 3 ^ 2 decreases to 
uj r = uj. At this resonant point the beam mode joins on to the L-O mode Ql7j], Beyond this point amplification ceases and the 
wave escapes as in the L-O mode. 

Let us estimate the gain factor using the approximation for T n for the cold plasma case, assuming 7 P and 7b do not 
change during the outflow, with uj r cx r^J 2 cx R~^l 2 The gain factor at a given uj is G = exp(2 f^ g T(uj, R)dR) = 
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exp[2(a;f ) oi?o/7b //2 ) :z; ( a; 2 ^ 3 — 1)]> where co p o is the plasma frequency at R = Rq, and x = lu /2lo p q^ 2 . The factor G is a 
maximum at x « 0.2 and, reverting to ordinary units, G max ~ exp(0.25w;,o-Ro/c7^ 2 )- For numerical estimates we use the 
following parameters: P = 1 s, B = 10 12 G, n p — Mncj = 10 2 x 6 • 10 10 cm -3 (near the pulsar surface), 7 P = 10 [l6| ], 
7b = 10 3 , rib/rip ~ 1, and Rq — 10R P JTT|]. This immediately gives that the maximum is achieved at u>/2ir « 500 MHz. The 
corresponding gain factor exp(G) ~ exp(30) « 10 12 , which implies efficient growth. Lower multiplicities would result in 
lower maximum gain frequencies. Lower 7;,, on the other hand, would make the instability more efficient. Better knowledge 
of the pulsar plasma parameters is required for direct comparison of our predictions with observations. The above estimate is 
invalid for low frequencies, corresponding to wave excitation at very large radii R > 10 2 i?o- At such radii one has uj r ~ fi/7 P 
where Q is the gyrofrequency, and then the infinite magnetic field approximation is no longer valid. Thus there should be 
a significant decrease of efficiency of wave growth at low frequencies. Specifically, using parameters chosen by Ref. || the 
spectrum should be cut off for u/u) r < 10~ 3 . 

Let us summarize the predictions of the proposed model. Quasi-transverse waves are generated efficiently in a wide frequency 
range below oj p qj^ 2 for small angles of propagation. The power spectrum of the escaping radiation is assumed cx G 2 , where G 
is the gain factor, and there is a maximum in G 2 at the frequency ui max ~ 0.1w p o7p^ 2 - The efficiency of wave growth decreases 
at both u> < w max and u> > w max , and the rate of decrease (steepness of the curve G(ui)) increases with the distance from w max . 
The implied form of spectrum is consistent with observations (e.g., Ref. [j]). 

We conclude that the nonresonant beam instability efficiently generates quasi-transverse waves in the radio range in a one-step 
process. The growing waves are in a beam mode well below the resonant frequency. As these waves propagate outward through 
the magnetosphere they evolve into the L-O mode when their frequency matches the resonant frequency. It is also the place 
where the growth ceases and, therefore, the spectrum forms. This emission mechanism plausibly reproduces the basic features 
of the observed radio spectra of pulsar, notably the existence of a maximum frequency with the spectrum steepening towards 
both higher and lower frequencies. However, a detailed interpretation of the observed power spectra requires that the emission 
mechanism proposed here be complemented with a statistical model for the emission. The observed spectra are obtained by 
integration over many pulses, and each pulse probably involves emission from a statistically large number of individual beams. 
Moreover, there remains a potential problem in explaining the lowest frequencies inferred for some pulsar emission: although 
the nonresonant model allows emission considerably below the resonant frequency, a more detailed investigation is required to 
determine whether the frequency range can extend low enough. A more detailed model for pulsar radio emission based on the 
mechanism proposed here will be presented elsewhere. 
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